Classical analysis of Bianchi types I and II in Ashtekar variables 



(N 



Gabriela Gonzalez * 
Department of Physics, Syracuse University, 
201 Physics Bldg., Syracuse NY 13244-1130 

Ranjeet S. Tate^ 
Department of Physics and Astronomy, University of Pittsburgh, 
100 Allen Hall, 3940 O'Hara St., Pittsburgh PA 15260 



gr-qc/9412015, UCSBTH-94-44 



We solve the complex Einstein equations for Bianchi I and II models formulated in the Ashtekar 
variables. We then solve the reality conditions to obtain a parametrization of the space of Lorentzian 
solutions in terms of real canonically conjugate variables. In the Ashtekar variables, the dynamics 
■ of the universe point particle is governed by only a curved supermetric - there is no potential term. 

' In the usual metric formulation the particle bounces off a potential wall in flat superspace. We 

^— ( . consider possible characterizations of this "bounce" in the potential-free Ashtekar variables. 

o : 



I. INTRODUCTION 



The Ashtekar variables for general relativity ^ are complex canonical coordinates on the real phase space, in terms 
I \ of which Einstein's equations become low-order polynomials. This has led to several simplifications and progress on 
J>. I various fronts ||^,^ . Classically, the use of the new variables has led to a greater understanding of the space of selfdual 
solutions of Einstein's equations |^ , as well as of solutions corresponding to degenerate metrics Q and some efforts 
I towards solving the classical constraints 1^ . In the canonical quantum theory, the use of new variables has led to the 
' discovery of large classes of solutions to the quantum scalar constraint equations and quantum states approximating 
given 3-geometries. Much of this progress has resulted from the simplified form of the constraints, in particular, 
the scalar constraint is a low-order polynomial, homogeneous and quadratic in the canonical momenta, albeit the 
Qv^ supermetric it defines is in general complex and non-flat. Thus, dynamical trajectories correspond to null geodesies of 
this supermetric and this simplification may allow one to solve the classical complex equations of motion, or at least 
to delve deeper into the structure of the space of solutions. 

This progress, however, has come at certain expense. In order to recover real Lorentzian general relativity, the 
^J^' new variables have to satisfy certain reality conditions. The configuration variable is a complex SU{2) connection 
. . : on the spatial 3-manifold S and the canonically conjugate momentum is a complex triad (of density weight +1). To 
recover real Lorentzian 4-geometries, one requires that the triads be real and that the real part of the complex SU (2) 
connection be the spin connection compatible with the triad. 

Can these reality conditions be solved? If one reintroduces the geometrodynamical variables one obtains the usual 
d formulation which is manifestly real, but in which the constraints and equations of motion are of a complicated 
form. The key question is: Can one first exploit the simple form of the constraints to solve the complex equations of 
motion and then attempt to solve the reality conditions? How does the space of Lorentzian solutions "sit" inside the 
space of complex solutions? Is there a simple characterization of the Lorentzian solutions in terms of real canonical 
variables? Clearly, such a real parametrization of the space of Lorentzian solutions can be obtained by starting from 
the geometrodynamical variables, however, our interest is to learn about the new variables themselves and the role of 
the reality conditions, since they appear to have greater potential for the full theory itself. 

This issue - that of the "reality" structure of the space of solutions- is also likely to be important from the point 
of view of the quantum theory. A criterion for the selection of a physical inner product on the space of solutions to 
the quantum constraint equation is to require that the reality conditions on physical observables be represented by 
Hermiticity conditions on the corresponding operators 0. Functions on the space of solutions are classical physical 
observables and thus may provide guides to constructing quantum Dirac observables. 

In the metric variables, the dynamics is that of a "point particle" (corresponding to a 3-geometry) moving under 
the influence of a potential in a real flat background superspace. In the new variables the dynamics of the system is 
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governed entirely by the curved supermetric, and since there is no potential term, the geometric structure of superspace 
assumes greater importance, and one would like to compare the behaviour in terms of connection dynamics to that in 
terms of geometrodynamical variables. Consider the Bianchi models which are certain homogeneous cosmological 
models of general relativity, with a finite number of degrees of freedom. In the geometrodynamical variables, the 
superspace can be taken to be flat and the particle bounces off potential walls which take different forms in the 
different Bianchi models. In Bianchi I the dynamics is that of a massless relativistic particle. In the simplest 
nontrivial case, Bianchi II, each dynamical trajectory bounces once off a single exponential wall and asymptotically 
the solution behaves like different Bianchi I solutions. In Bianchi type IX the potential has a more complicated shape: 
it is an expanding triangular well which confines the particle; and while the system has been shown analytically 
to be nonintegrable it is known that the particle bounces an infinite number of times off these walls The 
dynamics is a sequence of Bianchi I solutions -during which the particle is largely free of the influence of the potential- 
connected by Bianchi II solutions during which the particle bounces off one or the other walls. These bounces and 
their frequency play an important role in the numerical study of the dynamics and the presence of chaos in the system 
p| . In terms of the new variables, even though there is no potential, we expect that Bianchi II solutions behave 
asymptotically like different Bianchi I solutions, mediated by a transition corresponding to the bounce; and that the 
Bianchi IX evolution corresponds to an alternating sequence of Bianchi I and II solutions. Since in Bianchi II there is 
a single such transition between different Bianchi I solutions, it serves as a useful model to look for a characterization 
of this transition. Such a "unique" characterization would be useful for (numerical) studies of Bianchi IX in the new 
variables. 

In many ways Bianchi types I and II are prototypes for all the solvable homogeneous cosmologies, and in this paper 
we will focus our attention on the more complicated and interesting Bianchi type II. Both the models we will study in 
this paper have been solved in the metric variables by other means. We consider them here to learn something new 
not so much about the models themselves but rather about the utility of the new variables. 

The plan of the paper is as follows: In section 2 we review the canonical formulation of the Ashtekar variables in 
the context of the Bianchi cosmologies, following Ashtekar and PuUin |^ , and we discuss some structure common to 
all the Bianchi models. Section 3 is devoted to a study of the Bianchi type I model. While this model is quite trivial 
and has been quite thoroughly solved and understood in the new variables, the motivation here is to solve Bianchi I 
via a process which will be useful to solve Bianchi II, and to understand and familiarize ourselves with that process. 
In the first part of section 4 we will find the general solution to the classical complex Hamiltonian equations of motion 
for Bianchi II, and find that the space of solutions is parametrized by two canonically conjugate pairs of variables. In 
the second part of section 4 we will require that the triads be real throughout the evolution for all solutions. While 
the reality conditions are complicated in terms of the triad and the connection, they lead to simple conditions on the 
chosen reduced phase space coordinates: all the parameters turn out to be real, and thus the resulting reduced phase 
space description is simple. In section 5 we characterize the "bounce" during the transition from one asymptotic 
Bianchi I solution to another in terms of certain functions of the connection components. 



II. PRELIMINARIES 

Let us recall the Ashtekar variables for general relativity [^|J|]. These are the canonically conjugate pair {Ef, A^) of 
complex coordinates on the phase space of real general relativity, where Ef is the densitized triad and is a CS'C/(2) 
connection. These variables are related to the usual canonical metric variables via 

EtE"'^^'^" and Al^Tl-2Kl (1) 

where is the density weight +2 contravariant 3-metric on the spatial slices, FJj is the spin connection of the triad 
and is the (triad component of the) extrinsic curvature of the spatial manifold. 

Now let us consider the spatially homogeneous Bianchi models. The reduction of the Ashtekar variables to these 
cosmologies has been carried out previously ||^, p^. Here we will review the canonical formulation, following the 
conventions and notation of Ashtekar and PuUin]^] . The Bianchi models are spatially homogeneous models which 
admit a three dimensional isometry group which acts simply and transitively on the preferred homogeneous spatial 
slices. One introduces a basis (and its dual) of group invariant 1-forms lu^ , 1—1 — 3 which satisfy dui^ = —^C^jk^^'^ A 
Lu^ , where jk are the structure constants of the Lie algebra associated with the Bianchi group. Components of 
homogeneous tensors in this basis are spatial constants; e.g. the 4-metric in this basis is 

ds' = -N\t)de + Y.^u{t)u' <^oj' , 
IJ 
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where the lapse function N{t) and the 3-metric components qij are constant on the homogeneous spatial slices and 
are thus functions only of the time t. 

Class A models — distinguished by the vanishing of the trace of C — readily admit a Hamiltonian formulation . 
The structure constants for the Class A models can be written as 

jK — n'^^^ejjK (no sum over /), (2) 

where ejjK is the totally antisymmetric tensor and the set of numbers {n*-^^} — each of which can take the values 
0, ±1 — are used to classify the different models. In this paper we are interested in the Bianchi types I and II, which 
are both in Class A. 

For simplicity we will further confine our attention to the diagonal gauge, in which the off-diagonal components of 
the spatially homogeneous triad and connection vanish. The resulting phase space is 6 dimensional. In this diagonal 
gauge, the 6 canonically conjugate phase space variables are {Aj, —iEj}, with / = 1, 2, 3 and the symplectic form is 
given by 

-i^dEi AdAj. (3) 

For future reference note that the relationship between the 3-metric and the triads for the Bianchi cosmologies is 

E1E2E3 



qi 



E] 



(4) 



In the diagonal gauge the diffeomorphism and Gauss constraints vanish identically, and we are left with the scalar 
constraint: 



S = GNG''{A,n)pip, 



= -G 



N 



{AjAj - Sij{Aj)^) + ^n(^HeKij)^AK {-tEj){~tEj), 



(5) 



where Gn is Newton's constant, G^'^ is the (contravariant) supermetric on the configuration space of connections, 
and varies from one Bianchi model to another through its dependence on n^^\ Note that we have made a choice for 
the real lapse function: N ^ 1, thus N = {EiE2E-iY/'^ . 

In general relativity, classical dynamics is generated by the (vanishing) scalar constraint. Let A be an affine 
parameter along the orbits generated by the scalar constraint. Then the evolution of any function on phase space is 
given by / :— df /dX — {/, S\. Since the scalar constraint function in the new variables consists of only a kinetic term 
quadratic in the momenta, and is constrained to vanish, the dynamics it generates is that of a massless free particle 
in a curved spacetime with a (super)metric given by (||). Thus, in configuration space the dynamical trajectories are 
simply the null geodesies of the above supermetric. Even though the supermetric is complex, its geodesies and its 
null directions are both well defined. 

As we can see from the canonical transformation (0), the coordinates (Ei^Ai) are in general complex. In order 
to recover the real Lorentzian solutions, we have to impose "reality conditions" on the canonical variables. These 
conditions are: 

El 

1 d ( E^E2EzvS-^^\ (6) 



where the n^^^ defined in (|^) distinguish the various Class A Bianchi models. Both the scalar constraint (|5|) and the 
symplectic form (||) are real when evaluated on triads and connections satisfying the above reality conditions. Since 
the triads are also required to be real, the parameter along the real evolution generated by the (real) scalar constraint 
should also be real. We see that the reality conditions in the form i?/ e H, -E/ € IR — where the dot ■ indicates 
derivative with respect to a real parameter along the orbits generated by a scalar constraint smeared with a real lapse 
iV = 1 - are completely equivalent to the reality conditions in the "canonical" form (^. 

In the rest of this paper we will consider the Bianchi type I and II models, which are specified by rt*-^-* and 
n*^^^ — 1, v!-^^ — rS^'> — respectively. 
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III. BIANCHI I 



Now let us consider Bianchi type I, where the n'-^^^ = 0. The supermetric is then 



G'' = - 





A1A2 
A1A3 



A1A2 




AiA^ 





(7) 



In order to simphfy the calculations we will change coordinates on the space of connections, {A/} — s- {x^ — t, x, y} 
as follows: 



t 

X 

y 



A2A3 
GNA2A3/A1 

^2 Ma- 



in these coordinates the metric G is now diagonal: 

G'^ ^dtag{~t^,x^,y^), 

and the connection components are 

Ai = GNt/x 

A2 ^ Vty 
A3 - VtJy- 



(8) 



(9) 



(10) 



Note that there are many ways to diagonalize the supermetric, this particular choice of coordinates has been made 
in order to follow as closely as possible the change of coordinates which simplifies the calculations for the Bianchi II 
model. 

Since we are interested in the evolution of the triads, and will impose reality conditions on the triads, let us express 
them in terms of the new coordinates and momenta. Let pi = {pt,Px,Py} be the momenta canonically conjugate to 
the new coordinates t, x, y respectively, so that the triads are given by 



El 
E2 
E3 



XPx 
XPx 



ypy) 
ypy), 



and the symplectic form is 



= ^ dpi A dx' 



(11) 



(12) 



From (||) we see that under the transformation (^2,^3) (— A2,— ^3) the new coordinates are left unchanged, 



and there is the related ambiguity in the square roots in (10). We will resolve this ambiguity by taking the principal 
value of the square roots, and thus work with a quarter of the original phase space. The other parts are recovered 
by choosing the other signs of the square roots. We are justified in considering these "quadrants" separately since 
they are dynamically disconnected, i.e. the Hamiltonian vector field of the constraint is tangential to the boundary 
between any two of them. Furthermore, the spacetime geometries that result from the different choices of sign are in 
fact the same. 

Let us solve the (complex) equations of motion. For Bianchi type I the null geodesies are easy to find since the 
metric has three Killing vector fields: Ki — x''d/dx\ Let A € C be a parameter along the orbits generated by the 
scalar constraint, and let a;'(A) be an orbit in configuration space. Then, associated with the Killing vectors we have 
three conserved quantities : 



X'pi 

{tpt,xpx,ypy} 



G-jk m x>^ 
{-i/t,x/x,yly}, 



(13) 



We require the trajectories to be null, in order to satisfy the scalar constraint. This yields a condition on the above 
constants: 



- + + = 0. 



(14) 
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We can directly solve the equations of motion (|13|) for the orbits to obtain 

t{X) = toe-^^^ 

x{\) ^ e=-(^-^o) (15) 
y{\) = yoe'^y^, 

where the complex constants (to,Ao,?;o) correspond to the initial values of the coordinates {t,x,y) respectively (We 
have chosen Ao as a parameter instead of xq following the most convenient choice for Bianchi II) . 

If we substitute (15) into ( p^ ) and (pi]), we obtain the general complex solutions to the scalar constraint: we still have 



to impose the reality conditions. But before we proceed to the real solutions, we will first calculate the symplectic 
form on the complex reduced phase space F. We will do this in two steps. First, we evaluate il on the complex 
solutions (effectively taking as phase space coordinates the orbit parameters and the afhne parameter along the orbits 
{ct,Cx,Cy,to,yo,X}): 



^Ui(A) = J2idci Adlnx' 



(16) 



dct A dln<o + dcy A dlnyo + dX A (cfdct — Cydcy — Cxdcx). 

Next we pull-back to the constraint surface given by (|l4[), to obtain the symplectic form on the reduced phase space 

Q. = dct A dlnto + dcy A dlnyo. (17) 

Therefore we see that {ct, Cj,, Into, Inyo} is a set of (complex) canonical coordinates on the complex reduced phase 
space. We solve (^^ for the parameter as a function of the free parameters Ct,Cy: Cx = pr.val.^ — Cy. We will 

introduce in the solutions the discrete parameter e = ±1 to capture the ambiguity in sign of the solution to (p^). 

We now impose the reality conditions (^ on the complex solutions and find the phase space of the real degrees of 
freedom. Recall from the discussion in section 2 that this is equivalent to the requirement that the triads Ei be real 
for all "real" evolution, and that the parameter along the real trajectories should be real. Hence it is the real part of 
A which is the parameter along the real dynamical trajectories: if we write A = r + then t is the time, and (f) is 
a fixed real orbit parameter. The remaining reality conditions will lead in turn to conditions on the orbit parameters 
{ct, Cy, Into, Inyo, </'}• Let Aq = To + i(t>o, and identify tq as the real initial time. 

We define new orbit parameters 

To = to exp(-ict(?!)) and = yo exp(iCy(/)). (18) 
From the expression (O) and the solutions (113) we obtain the solutions for the triads 



cc cc 
El = -i-^ = ^ exp((ct + eCx)T - ec^To + iecx{(f> - 0o)) 

.Ct + £Cx + Cy .Ct + <^Cx + Cy 
E2 =1 ^=--^exp((ct -Cj,)r/2) 19 

Es. ^i- — ^ — ^ i{ct + ecx-cy)./^ expect + Cy)T/ 2). 

A3 V 

Let us impose the rest of the reality conditions. Since we only have exponential functions of r, we see that d In Ej / dr 
are real if and only if all c,; are real (and hence Cj > c^). Next, we see that E2 and Ej, arc real if and only if in addition 
To, Yo are real and satisfy Toyo < 0. Then, Ei will be real only if exjp{iecx{4> — (f'o)) — i*, from which we conclude 
that cos{cx{<j) — 4>o)) — 0. 

Collecting all reality conditions, we conclude that the triads are real if and only if 

Ct, Cj^ e IR and Cj > Cy 
ro,yoeIR and Toro<0 (20) 
and cos{cx{(f> — 0o)) = 0. 

Note that To corresponds to a choice of the initial value of x and has been left as a free parameter, but (f> and (j>o have 
disappeared as orbit parameters, fixed by the reality conditions. The real solution for the triads is then simply 

El = ± ^ exp((ecj; + ct)T - cCxTq) 

E, = ^^+!^^^ exp((c, - Cy)r/2) (21) 



E3 = {ct + ecx - Cy)\J exp((cf + Cy)T/2). 
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Note that the sign ambiguities in the above solution for the triads lead to distinct triad solutions. However, the 4- 
dimensional spacetime geometries (quadratic in the triads) are unaffected by the choice of sign and thus the solutions 
are physically equivalent. 

When the reality conditions are satisfied, the connections are pure imaginary: 

Ai = ±zGVToe^"-"°exp(-(ct + ec,)T) 

A2 = ^^/5^^exp(-(Q-Cy)T/2) (22) 
^3 = ^v/-^o/ro exp {-{ct + Cy)r/2) . 

All the connections have a simple exponential behavior, vanishing (at different rates) as c^r — > 00. Note that there 
are sign ambiguities in all three connection components arising due to either the phase ecx{(l)— (f>o) or the square root. 
These sign ambiguities also lead to different solutions which however correspond to the same spacetime geometry. 
Finally, the pull-back of the symplectic form to the space of real solutions still looks very much as in (pT|): 

Cl = dct A dlnTo + dcy A dlnYo. (23) 

Thus, the chosen coordinates on the reduced phase space are real as well as canonically conjugate. 

From the above solutions and (Q), one can obtain the components of the metric. The metric is the Kasner solution 
0: 

ds^ = -dt^ + t^P' dx^ + t^P^ dy^ + i^P^* dz^, (24) 
where we have rescaled the time function 

g(«.+2cOr ^ ^ (25) 

and the set of Kasner parameters p = {pi,P2,P3} are functions of the orbit parameters 

{ct,Cx, Cy} p ^{~ecx, ecx + Cy - Ct,eCx - Cy - Ct}/{ec^ + 2ct). (26) 

The properties of the Kasner parameters X) = = 1 are derived from (|l^) . Belinskii et al introduced a 

single parameter u G [0, 1) as a solution to the above constraints, where (1 + m + 'u^)p(u) — {— m, 1 + m,w(1 + w)}. 
In terms of the parameters we have used, v? = {\ct\ — \cy\)/{\ct\ + |cj^|) for the ratio \cy/ct\ in the range [0, 1]. The 
solutions with parameters |ct/cj,| — 1 correspond to flat Rindler spacetime (i.e., Minkowski spacetime in an accelerated 
frame). The "trivial" solution with all Ci = is, of course, Minkowski spacetime. In the non-trivial cases, the initial 
singularity is approached in the limit — CfT — s- 00. 

IV. BIANCHI II 

In this section we will solve the classical dynamics of the Bianchi type II model, which is specified by the structure 
constants n^^^ = 1 = C^23- We will repeat the procedure step-by-step as for Bianchi type I. 
The Bianchi type II supermetric is 



A1A2 AiA^ 

A1A2 -I- Ai/Gat 

A1A3 A2A:i + Ai/GN 



(27) 



As for Bianchi type I, we will change coordinates to a set adapted to the Killing symmetries of the supermetric, 
{Al} — + {x^ — t, X, y} as follows: 

X = GNA2A3/A1 (28) 
y = A2/A3. 

In these coordinates the metric G is now diagonal: 

G'^ =dtagi ^.t^,x{x~l),^-y'). (29) 

X ~ 1 X 

The connection components are expressed in terms of the Killing coordinates by 
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= G Nt/{x-\) 

A2 = ^txy/{x ~ 1) (30) 

A3 = VW(y(a^-i))- 

(The ambiguity in the sign of A-^ we resolve as for Bianchi type I; see the discussion after (|l0|).) Let "Pi — {vt^Vx^Vy] 
be the momenta canonicaUy conjugate to the new coordinates x, y respectively. Then the triads are given by 



2 



A2 

i 



(^^— Y*Pt + ^Px + ypy) (31) 



E3 = -r-i Ttpt + - ypy). 

Ji-l X — L 



i3 

For Bianchi type II the supermetric has only two linearly independent Killing vector fields. However, since the 
space is 3-dimensional, the system is separable, and we will be left with three ordinary differential equations to solve. 
The two Killing vector fields are: Kt = td/dt and Ky = yd/dy. As before, let A £ C be a parameter along the orbits 
generated by the scalar constraint, and let a;*(A) be an orbit in configuration space. Then, associated with the Killing 
vectors we have two conserved quantities : 

1 — X i 

ct = tpt = tGtt t = - 

- I (32) 

Cy ^VPy =yGyyy^—— - 
X -f- 1 y 

Let us define Cx similar to the Bianchi I parameter: 

Cx =pr.val.^cj - cl- (33) 

The constraint (i.e. that the trajectories be null) yields an equation for x 

x'=x'cl + cl. (34) 

This equation can be readily solved for x{X): 

x(A) =e^sinh(c,(A- Ao)), (35) 

where e = ±1 corresponds to the two possible choices of sign in the square root of (|^. 

Substituting the above solution a;(A) in the equations ( |32| ) for the conserved quantities we obtain ordinary linear 
differential equations for t{X) and y(A) which can be solved to yield the complete equations of motion: 

t{X) = toe-^^^^-^t^i^ 

Cy X Cf 

^ _c*ACa:COsh(c:r(A - Ao)) + Ct sinh(c2:(A - Ao)) 

= n~~r\ — (36) 

Cy cosh [Cx (A - Ao) ) - act 
x^-Cy _ c„A'^Qsh(c^(A-Ao))-£ 



y{\) = y^e^y^ ^ = yoC 



sinh {cx{\ - Ao)) 



where the complex constants ioj2/o correspond to the initial values of the coordinates t^y respectively. 
A lengthy calculation shows that the symplectic structure evaluated on the solutions is simply 

f2 = dct A Into + dcy A Inj/o, (37) 

from which it is clear that {ct, Cy, Into, Inyo} is a set of canonical coordinates on the complex reduced phase space. 
The calculation of the symplectic form on the reduced phase space is considerably simplified when we note that 
the pull-back of the symplectic structure to the constraint surface has two properties:zj it is degenerate -with the 
degenerate direction along the Hamiltonian vector field of the constraint- and ii) its Lie derivative along the above 
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vector field vanishes. Hence, there can be no terms in Cl containing dX, and the coefficients of the remaining terms 
will be independent of A. Thus, we can simply evaluate the symplectic structure on some A = const cross-section of 
the constraint surface. 

We now impose the reality conditions (|^) on the complex solutions and find the reduced phase space of the real 
degrees of freedom. As before we can identify r = i?e(A) with the real time, and require that the triads and their 
derivatives with respect to t should be real. It is more convenient to (partially) solve the reality conditions in terms 
of E2 ■ E3 and E2/E3. First, we substitute for the momenta pi in terms of the velocities in the expressions (|3l|). 
Various simplifications then lead to the following expressions: 

Erir) E2/E3 = ^ ^^-^ 



Ep{T) := E2- E3 = -—c^cye^'-, 
-to 

where we have defined the new orbit parameters To, lo exactly as for Bianchi I (^8|). Now d\n{Er)) / dr and dln(E,p))/dT 
are real if and only if 

Cy^cte-R, (39) 

and the reality of Ep, Er is then equivalent to 

cjYo,ajToeWi, (40) 

where Cx is either real or imaginary. In order to guarantee that £'2,i?3 are themselves real, and not possibly pure 
imaginary, we also require that Er ■ Ep > 0, this condition on the orbit parameters is 

Cy{ct+Cy)YoTo<0. (41) 

Now let us consider the reality of Ei . With the above reality conditions on E2 , E^ imposed, the solution for Ei is 

— ^ 

Eiir) = T^e"''^ (ct sinh {cx{t - tq + - 0o)) - cosh {cx{t - to + i(0 - 0o))) , (42) 
Jo 

where, as for Bianchi type I we have defined Aq = tq + i4>Q. We will impose the reality conditions on Ei,Ei at r = tq. 
Using the previously solved reality conditions on the orbit parameters, the remaining reality conditions are reduced 
to 

£;i(To)eIR-^ —su\{cx{<j3-(t)Q)) + icos{cx{<j3-(t)Q))e'B. 

Cx (43) 
and i?i(ro) e IR -^=^ Cx sva{cx{(t) - 4>o)) + ict cos{cx{(l) ~ M) ^ M. 

Now we know that Cx is either real or pure imaginary; in either case the first terms in the above equations are real, 
and the second terms are pure imaginary and hence should vanish. Thus, we should have cos{cx{4> ~ 4'o)) = 0; which 
has no solutions for imaginary c^, and we conclude that c^: € IR and thus Cj > c^. Putting this back into |4^, we see 
that these (necessary) conditions are sufficient to guarantee that Ei is real for all r. 

Collecting all the reality conditions, we conclude that the solutions for the triads in Bianchi type II will be real if 
and only if 

Ct, Cj^ e IR and > Cy 
ro,roeIR and CtCyToYo<0 (44) 
and cos{cx{4> — 0o))) = 0. 

(In the inequality (|4l| ) the sign of (ct + Cy) is determined by the sign of ct. Note also that tq, which corresponds to 
the initial value of x, is left undetermined by any of the reality conditions, and is a free parameter.) 
The triads evaluated on the real trajectories are: 

El = icfrb *5xp(ctT) cosh(c;j;(T - To) - 0) 

■E^2 = \/ T^,Cy{ct + Cy) exp(ct - Cy)T/2 (45) 
E3 = ^ =^Cy{ct - Cy) exp(ct + Cj,)r/2, 
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where (j) =arctanh(ca;/ct). 

If we evaluate the metric, we find the solution quoted by Belinskii et al [p^ (as an approximate solution for a 
"Kasner epoch" transition in Bianchi IX) with the same identification of parameters we made in Bianchi I in (|2^), 
with e = +1. It is also, of course, the solution first quoted for Bianchi II models by Taub [|5|, with a similar 
identification of parameters k = Cx,ci — ct — Cy and C2 — Ct + Cy . 

The connections Aj evaluated on the real trajectories are : 

Al = -GN{cx/cy)Toexp{-ctT){l + ief{T)y^ 

Al ^ Toroexp(-(Q - c,Jr)(l - iesinh(c,(T - tq))) (1 - ze/M)"' (46) 
Al - -(To/Fo) exp(-(Q + Cy)T)(l + iesinh(c,(T-ro)))(l-ie/(T))-\ 

where /(r) — 5 smh{cx{T — tq) — (p), with 5 = sgn{ctCy) . 
As T — !■ ±oo, the connections behave respectively as: 

Al -ieGNTaeyiY>{±CxTa)/{lTctlcx) exp(-ct ± c^r) 



A2 ^ i^/WoYo\\{ct±Cx)/ Cy\ eXp(— (Ct — Cy )r/2) (47) 

Cy\ eXp(— (Ct + Cy)T/2). 

We can see that these limits correspond to Bianchi I solutions, that is, Bianchi II solutions approximate (different) 
Bianchi I solutions as t ^ ±cx). These Bianchi I solutions differ in the value of the discrete parameter e and 
renormalizations of Tqi 

As in the case of Bianchi type I, the pull-back of the symplectic form to the space of real solutions is: 

(l = dctAd\nTfi + dcyAd\nYQ, (48) 
where Ct, Cy,To,Yo are Dirac observables for the theory. 

V. CHARACTERIZATION OF THE "BOUNCE" 

When solving Einstein's equations in the Hamiltonian formulation for the "diagonal models" , the usual geometro- 
dynamical description ||l^] is in terms of variables {(3"^, pa}, with the index A taking values 0, ±1 where is the log of 
the 3-volume and are the anisotropics. The Hamiltonian constraint has a term quadratic in momenta {rj^^pAPs) 
and a potential term which for all Bianchi types (except Bianchi I) has exponential dependence on the 2-dimensional 
plane. Therefore the dynamics is that of a particle (the "universe point" ) moving on a flat background in the given 
potential. In the regions where the potential vanishes the particle moves essentially free along a (null) straight line as 
if it were a solution to Bianchi I equations (Kasner epoch), until it encounters a potential wall and bounces back along 
a (reflected) straight line. In Bianchi II there is only one "wall" in the potential, and therefore the world-particle 
bounces once between the initial singularity and the expanding evolution. In Bianchi IX models, there are 3 walls 
and in general the particle keeps bouncing within the walls, approaching the singularity as i — s- 0. 

This classical behaviour was flrst analyzed by Belinskii et al |l^]. They use the Bianchi II solution as a transition 
between Kasner epochs in the general Bianchi IX solution. In each Kasner epoch, the solution is approximately the 
Kasner solution, where the metric component in one spatial direction decreases with time (it has a negative Kasner 
exponent) and the other two increase with time (they have positive exponents). In the transition between Kasner 
epochs, the negative power of time is transferred from one spatial direction to another, so one metric component reaches 
a minimum and another a maximum, while the third one increases monotonically during the transition. Bianchi II 
models have a single transition, and the "bounce" is then deflned as the point along the trajectory when one of the 
metric components reaches its maximum (which is different from the point when the other metric component reaches 
its minimum). 

Now, if we consider the Hamiltonian formulation in terms of the new variables, the Hamiltonian constraint contains 
just a term quadratic in momenta. Therefore, the dynamics can be interpreted as that of a free particle moving 
on a null geodesic in the complex 3-dimensional space {^7}, where the supermetric is given by G^^ in (^). Since 
there is no potential, there is neither a wall nor any particular coordinate time to identify with the "bounce" time 
that characterizes Bianchi II in the geometrodynamical formulation. If we translate the bounce time deflned in the 
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geometrodynamical variables, we find that it is the point on the trajectory where the function t(r) vanishes, that is, 
when the KiUing vector Kt vanishes. However, since this description is in terms of KiUing symmetries that do not 
exist in some of the higher Bianchi models, this characterization of "bounce" times is not useful for the study of the 
other models. 

Since the dynamics is determined completely by the supermetric G^'\ one could look for an invariant description 
of the bounce in terms of the geometry of the curved superspace in which the "universe point" moves. Any scalar 
quantity that we construct from the supermetric is a function only of the coordinate x(r), which when evaluated along 
real trajectories depends only on the function cosh(ca;(T — tq)). Therefore there is a symmetry that one can identify: 
(r — To) — > — (r — To). Again, however, this symmetry does not survive when we consider other Bianchi models. 
Hence, we would like to find a way of describing a turn-around point that does not depend on the use of symmetries 
particular to the Bianchi H solutions. One possibility is to look for features such as maxima or singularities in (the 
absolute values of) various superspace curvature scalars like the Ricci scalar RjjG^"^ , or RijR^"^ . Unfortunately, all 
the real metric scalars up to second order in the Ricci tensor have several extrema as functions of t when evaluated 
on the real trajectories, and therefore do not uniquely identify the bounces in higher Bianchi models, like type IX. 

There is however, another option. The Bianchi H solutions behave asymptotically (r —^ ±00) like specific Bianchi 
I solutions. On real Lorentzian Bianchi I solutions, the connection components are pure imaginary, exponential 
functions of r. Hence we can look for a characterization of the "transition region" , by studying the ratio of real to 
imaginary parts of the solutions for the connection components in Bianchi type H: we know that this ratio tends to 
zero in the asymptotic regions as r ^ ±cxd, and since the Bianchi H connections are in general complex, this ratio is 
non-zero at finite r. The functions $/ — Re{Ai) / Im(Ai) turn out to be very simple: 

= eS sinh(ca;(T - tq) - <j)) 



-2 



ecosh(c2;(r - To) - 0/2)/sinh((/)/2) if 5 = +1 

-esinh(c:r(T-To)-</)/2)/cosh(</)/2) if 5 = -1 (49) 



$3' = esinh(c:r(T- To) - (/)/2)/cosh(</>/2) if 5^+1 
= -ecosh(c^(T- To) - 0/2)/sinh(0/2) if^ = -1, 

where we recall that S — sgn{ctCy) — —sgn{ToYo) and (f> — arctanh(ca;/cf ). Each of these functions has a unique and 
very obvious "transition point" (where they either have a maximum or diverge) at either Cx{t — tq) = (j) (for $1) or 
Cxir — To) = 4>/2 (for $2, $3). The zero of is exactly the "bounce time" as defined in [lOj. Thus, the "Kasner 
epochs" in the Ashtekar variables are identified with pure imaginary, exponential connections and the "bounces" are 
either maxima or divergences of the (tangent of the) phase of the connection components, or in other words, the most 
extreme departure from Bianchi I connections. Both before and after the bounce, the connections have an asymptotic 
exponential behavior, but with different logarithmic velocities. This description is very similar to that by Belinskii 
et al in terms of the behaviour of metric components. As in the geometrodynamical description where the maximum 
of one metric component does not coincide with the minimum of the other (forbidding an invariant definition of the 
bounce time), here too the divergences and maxima in the phase of the connection components do not coincide. 



VI. CONCLUSION 



Let us briefiy recapitulate the results we have obtained. For both Bianchi types I and H, we first solved the complex 
equations of motion by finding the null geodesies of the complex supermetric on the space of connections. The space 
of solutions is parametrized by 4 complex variables, which form 2 canonically conjugate pairs. Next, to find the real 
Lorentzian solutions, we required that the triads are real throughout the real evolution. This condition is satisfied 
when the parameters on the space of solutions are all real and satisfy in addition certain nonholonomic constraints. 
(Note that we have made a simplifying choice of variables in order to obtain a real canonical parametrization of the 
reduced phase space.) 

All Bianchi H solutions approach a (different) Bianchi I solution asymptotically in the past or the future, and 
deviations from Bianchi I like behaviour and the transition from one asymptotic Bianchi solution to another occurs 
at some finite time. Now for Lorentzian Bianchi I, Re{Ai{t)) = 0, whereas the Lorentzian connections for Bianchi 
type H are in general complex. So it is particularly nice to characterize the bounce in Bianchi H as a deviation from 
Bianchi I like behaviour via maxima or divergences in <&/ = Re{Ai{t))/ Im{Ai{t)). We have seen that for any given 
solution, one of the phases $7 has a maximum while the other two have a singularity. This does uniquely characterize 
the bounce and could be used for the counting of bounces in the numerical study of Bianchi IX in the new variables. 
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Recall that we failed to find a description independent of the coordinates on the (minisuper)space of connections. 
In retrospect it is not surprising that there is no characterization in terms of some supermetric curvature scalars, since 
the Lorentzian solutions have not been described as null geodesies in some real supergeometry, though the complex 
solutions are null geodesies in a complex superspace. 

One can proceed to construct the reduced space quantum theory for these models. This construction is straight- 
forward since the reduced phase space is coordinatized by two pairs of real canonically conjugate variables. However, 
from the point of view of learning something about the Dirac quantum theory in terms of the connections themselves, 
this is not a very useful approach. In particular, one does not see clearly the role of the complex reality conditions in 
finding an inner product on the space of solutions to the quantum scalar constraint equation ||^ , or whether this is 



even feasible. These issues will be explored in a later paper |17|, in which we construct the complete Dirac quantum 
theory for the Bianchi II model. The classical Dirac observables (cf , Cy, Tq, Iq) we have constructed here do have 
quantum analogs, and the Hermiticity conditions on them serve to select an inner product uniquely on the space of 
solutions to the quantum scalar constraint equation. 

VII. APPENDIX: SPECIAL CASE SOLUTIONS 

For special values of the parameters Cy,ct, the limits of the general solutions we found in section 4 are not themselves 
solutions to the limiting equations of motion. In order to find the correct solutions at these limiting values, we have 
to solve directly the reduced equations of motion obtained by taking the appropriate limits of the general ones. There 
are three special cases and we consider them one by one. 

Case 1: Ct — 

In this limit, the equations of motion (|3^ , |3^ ) reduce to 

i =0 

y/y =Cy{l + l/x) (50) 
i:2 =c2(l-a;2). 

The solution to the above equations is 

t{\) = ietci 

/^\ ■ c X sin(cj/(A - Ao)) 

yiX) ^lyoe^ — tt — -^—tt (51) 

' ^ l + ecos(c.y(A- Ao)) ^ ' 

x{X) = esin(cj^(A — Ao)). 

Note that the order and the form of the limiting equations are the same as in the general case, and thus we expect 
that the above solution to the limiting equations can be obtained by taking the limit (ct '—^ 0) of the general solution 
( ^ , ^6| ), as can be confirmed by direct calculation. Thus we can directly take the limits of the reality conditions (|44|), 
and conclude that Cy = 0. The triads Ej^s in this case (with ct — Cy — 0) vanish, and the connections A/'s are 
constants. 

Case 2: Cy = 



The equations of motion in this case are: 



i/t ^ctx/{l~x) 
y =0 (52) 



• 9 2 2 

X = CfX . 



and the solution is 



a:(A) =e'='(^-^«) 

y(A) =yo (53) 
t{X) ^ to{x{X) - 
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If we take e = —1, then E2 = = and Ei = ict/{GNto), and the connections are Ai = to, A2 = yo, A3 = 
y/tauo exp(— ct(A — Ao)/2). The reahty conditions only require that Ct/io is pure imaginary. 
If we take e = +1, then the triads are: 

E, ^ ^£^(i_e2-*(A-Ao)) 

2fi" (54) 
y/toyo 

The reahty conditions then require q, Uq, j/q G IR , cos(cf ((/> — (/jq)) = 0, and ito2/o sin(ct(^ — (po)) > 0. Let us define 
d = sin(cf (0 — 0o)) = ±1, and Tq — iSto, Yq — 7/0 (so TqYq > 0). The Lorentzian solutions for triads and connections 
are, then: 

E, = _^(i-e2^'(--o)) 
9° 



E2 = Es/Y, 







A V A ■ l^rTT exp(ct(r-ro)) 

A2^YoA-i = ly/ToYo- --. 

1 - idexp[ct{T -To)) 

Notice that although the triads have the same form as (^) (ignoring constant multiplicative factors, and taking 
Cy = 0,Cx — ct in the exponentials), the connections do not have the same behavior as in (|4^). In other words, the 
limit Cy —f can be taken smoothly in the metric variables, but not so in the connection space. 



We set Cy — dct, where d = ±1. The equations of motion are: 

t/t — ctxj (1 — x) 

y/y =5ct{l + l/x) (56) 



and yield the solution 



x(A) =ect(A-Ao) 

y{X) ^yoiX-XoY'e'^^^ (57) 
t{\) -to(l-eQ(A- Ao))-'e-^*\ 

We see that S —^ —S is equivalent to y y^/y, and this in turn is equivalent to {A2 < — > 2/o^3i ^2 < — * yo-£'3}i so we 
only need to consider in detail the case 5—1. 
If e = +1, the triads are: 

E, = ^e^^\l-c^{\-\o)) 
_ pri + Q(A-Ao) 

E, = y^e^'Ml + Q(A-Ao)). 

We can see by inspection that the only solution to the reality conditions is the trivial one, ct = 0, (vanishing triads 
and constant connections). 



If e = — 1, then the triads are: 



El = Q 

. ~cr /l + ct(A- Ao) 
E2 = I 



toyo\ 1 - Ct (A -Ao) (59) 
p ^ . fc^j^ / l + ct(A-Ao)" 
' to A-AoV l-ct(A-Ao)- 
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Here, the reality conditions are satisfied when (j) = = and Ct-,to, and yo are real and satisfy Cftoyo < 0. 
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